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Abstract. We study charged spin-1/2 particles in two dimensions, subject to a 
perpendicular non-Abelian magnetic field. Specializing to a choice of vector potential 
that is spatially constant but non-Abelian, wc investigate the Landaii level spectrum 
in planar and spherical geometry, paying particular attention to the role of the total 
angular momentum J = L + S. After this we show that the adiabatic insertion of 
non-Abelian flux in a spin-polarized quantum Hall state leads to the formation of 
charged spin-textures, which in the simplest cases can be identified with quantum Hall 
Skyrmions. 
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1. Introduction 

It has long been known that the motion of charged particles in a plane and in the presence 
of a perpendicular magnetic field is highly special. In classical mechanics particles 
exhibit cyclotron motion, while quantum mechanics leads to the Landau level band 
structure [Ij , with each level providing a macroscopic number of one-particle states that 
are strictly degenerate in energy. In this setting free electrons can form integer quantum 
Hall states, with intricate topological properties resulting in a quantization of the Hall 
conductance. Adding interactions leads to a plethora of remarkable many-body states, 
collectively known as fractional quantum Hall states. Such states are known to admit 
excitations with fractional charge and fractional statistics. The possibility that specific 
quantum Hall states have excitations with non-Abelian braid statistics has opened the 
exciting perspective of applications in the realm of topological quantum computation 

A setting that is physically very different but mathematically similar to that of 2D 
electrons in a perpendicular field is that of rapidly rotating cold atoms. In this analogy, 
the vorticity of a rotating liquid is akin to magnetic flux in the electron system. The 
limit of rapid rotation leads to a Landau level structure and one expects the formation 
of atomic quantum Hall states (see |3] for a review). The road towards experimental 
realization of atomic quantum Hall states is extremely challenging, but the lowest 
Landau level (LLL) has been reached ^ and there are indications of the formation 
of incompressible quantum liquids in small clusters [5]. 

The cold atom setting allows for yet another variation on the same theme: cold 
atomic gases subject to external time- dependent potentials that are such that they 
mimic the effects of magnetic fields or, equivalently, of rotation. Recent proposals for 
implementing such artificial gauge fields can be found in [U [TJ [H [U [TOl [TTl [12]; a 
successful experimental realization was reported in [13] . 

Very interestingly, the latter setup is flexible enough to allow for a generalization 
of regular, Abelian, magnetic fields, corresponding to the U{1) gauge symmetry of 
Maxwell's theory, to gauge fields pertaining to non-Abelian symmetries. The simplest 
version of this idea is the case of spin-1/2 particles subject to an external gauge potential 
for a non-Abelian U{2) symmetry [111 [151 [161 [17] , see [18] for the case of spin-1 particles. 
Remarkably, very similar if not identical Hamiltonians arise in the original setting of a 
2D electron gas in a perpendicular magnetic field, if Rashba [IH] and Dresselhaus [2U] 
spin-orbit coupling terms are taken into account. This connection can be exploited to 
study the physics of the quantum spin Hall effect in the context of cold atoms [211 [22] • 

With this, there is ample reason to study non-Abelian external gauge potentials 
with special emphasis on what is new as compared to the Abelian case. For particles 
confined to a lattice geometry, the non-Abelian case leads to an interesting generalization 
of the 'Hofstadter butterfly' fractal 1-particle spectrum, dubbed 'Hofstadter moth' [15] . 
The Landau level problem and the possibility of realizing (fractional) quantum Hall 
states with particles subject to non-Abelian gauge fields has already attracted a lot of 
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attention [Ml [13 |23l El] . 

In this paper we analyze non-interacting spin- 1/2 particles in an external non- 
Abelian field and reflect on some of the fundamental differences with the Abelian setup. 
We analyze and solve the Landau level problem in spherical geometry, highlighting the 
fundamental role of the total angular momentum J = L + S, which commutes with the 
Hamiltonian. In this setup, the non-Abelian field penetrating the sphere agrees with 
the asymptotic (large radius) limit of the non-Abelian magnetic monopoles first written 
by 't Hooft and Polyakov [25j. One reason to focus on spherical geometry is that this 
geometry is known to be particularly useful for the purpose of a numerical study of 
many-body states arising upon adding interactions to the Landau level problem 

We then proceed to the process where a non-Abelian flux is inserted in a background 
of otherwise Abelian flux. In our simplest case, = + Sz remains a good quantum 
number during the flux insertion. This suggests a prominent role for transitions where 
particles flip their spin while at the same time changing their by (plus or minus) 
one unit, meaning that they jump to an adjacent Landau level orbital of the Abelian 
problem. This is indeed what happens: starting from a spin-polarized integer quantum 
Hall state, inserting non-Abelian flux at the origin leads to a state where, depending 
on their distance from the origin, particles have an amplitude for changing their spin 
and moving radially in or out. The resulting state is a spin-texture of unit electric 
and topological charge, which is easily identified with a quantum Hall Skyrmion. More 
general external fields lead to more intricate textures. 

In the literature, thought experiments involving insertion of (Abelian) flux are often 
invoked as a probe of the characteristics of the quantum phase of a many body system. 
One spectacular example is the argument by Laughlin that insertion of a unit flux 
through a gapped medium with Hall conductance an = ve^ /h leads to the nucleation 
of an excitation with fractional charge e* = ve [2Z]. In the case of the quantum spin 
Hall state, insertion of flux leads to spin-full excitations at the sample edges [28]. One 
motivation for the present study has been the wish to extend these considerations to 
the case of non-Abelian flux. 

We briefly comment on possible experimental realization of our ideas. In the case 
of electrons, the integer quantum Hall reference state is readily available, and the non- 
Abelian flux can in principle, in the presence of spin-orbit coupling, be generated via a 
perpendicular electric field (see section 3 below). Integer quantum Hall states for cold 
atomic fermions are expected to arise in rapidly rotating systems |29j or in laser-induced 
artificial gauge fields [30] . The paper [15] has indicated how non-Abelian external gauge 
fields can be generated with the help of external lasers; if such laser pulses can spatially 
focussed, this will result in pulses of non-Abelian flux similar to the ones we propose 
here. Clearly these external fields will deviate from our expressions in their details, 
we expect however that the spin-textures that we predict here are relatively robust. 
Reading out spin-textures of effective spin- 1/2 atomic states is possible; see for example 
[21], where a Skyrmion lattice for effective spin-1/2 bosons was imaged. 

This paper is organized as follows. We start (section [2]) with a general exposition on 
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non-Abelian gauge fields. Next (section |3]) we discuss the non-Abelian Landau problem, 
with particular emphasis on spherical geometry. Section |4] discusses the response of spin- 
full particles to the insertion of non-Abelian flux in an otherwise Abelian background. 
We end with some conclusions in section [5j The appendices A-D present further details 
and background material. 

2. On non-Abelian gauge fields 

Before we start to investigate the Landau problem, we start by a quick review on 
non-Abelian gauge fields, highlighting the differences with the Abelian case. We are 
interested in the quantum problem of a particle coupled to an external non-Abelian 
gauge field A. In a cold atom experiment these fields are controlled by external lasers. 
Therefore in the present article we are not concerned with the dynamics of these gauge 
fields, which are treated as control parameters. The Hamiltonian for a non-relativistic 
particle of mass m in an external magnetic field is 

In the non-Abelian setup each component of the vector potential A is a matrix 

A = AxUx + AyUy + A^Uz = A^Ua . (2) 

The corresponding magnetic field (or curvature) B is 

B = V X A-lAx A , (3) 

h 

or equivalently in components: 

Ba = eabc{dbA^-^-Ai,A)j . (4) 

The first term is the usual curl, while the second term A x A vanishes identically for 
Abelian gauge fields. However it is non-zero in the generic non-Abelian the 
components Aa of the potential may not commute. In contrast to the Abelian case, 
even a uniform potential A can produce a non-zero magnetic field. Introducing the 
covariant derivative is D = V — |A = |(p — A), the magnetic field takes the following 
compact form: 

B = Dx A . (5) 

A gauge transformation is simply a unitary transformation U and a change of A, such 
that D transforms covariantly: 

A UAU^ + ihUVU^ . (6) 

Then the following quantities transform covariantly (note that the magnetic field is no 
longer gauge invariant): 

D UDU\ B UBU^ , (7) 
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and the Hamiltonian ([T]) is clearly a covariant quantity: 

H=—{p-A)^ = -—D\ (8) 

In the Abelian case, two field configurations A yielding the same magnetic field are 
necessarily equivalent up to a gauge transformation (on a simply connected space). 
This is no longer the case for non-Abelian gauge groups. In particular it is known from 
[32j that there are two gauge-inequivalent kinds of non-Abelian vector potentials which 
produce a uniform B field (although this statement is completely general, we illustrate 
this for B = 2(7 zUz)'- 

• a commuting field with linear potential A = ^B x r = a^i—y, x, 0), for which only 
V X y4 contributes to the field strength. This case is Abelian in nature, in the sense 
that all components of the potential vector commute with each other. 

• a uniform non-commuting potential, for instance A = {—ay,ax,0), such that 
B = —j^A X A. This is only possible in a non-Abelian gauge group. 

Although these two kinds of potential give rise to the same magnetic field, they lead 
to completely different physical properties. While the first kind induces a magnetic 
length, Aharanov-Bohm effect and Landau level discrete spectrum, the second one has 
no Aharanov-Bohm effect, and the spectrum of a particle in such an external field 
configuration is continuous. 



3. The non-Abelian Landau problem: particles in perpendicular uniform 
non-Abelian external field 

In the present article we focus on the simplest case of non-Abelian gauge fields, when 
Aa and Ba are 2x2 hermitian matrices. The gauge group is then U(2) = f/(l) x SU{2) 
and decomposes into: 

• an Abelian U{1) part, namely the fields proportional to the identity matrix I, 

• and a non-Abelian SU{2) component, whose fields are linear combinations of the 
Pauli matrices aa- 

The U{2) case is a natural choice as it is the simplest case allowing non-Abelian gauge 
fields. However there is a deeper reason to focus on U{2) gauge fields. Very similar 
physics can arise in a 2D electron gas when taking into account relativistic corrections 
in the Pauli-Schrodinger equation, such as the Thomas term 

«r--^,S.{E.f). (9) 

This spin-orbit coupling term plays a crucial role in spintronics ^33j, and it mimics the 
effect of a non-Abelian gauge potential: 

Ar^Exa. (10) 

In this section we study the quantum problem of a non-relativistic particle confined to a 
two dimensional manifold in the background of a uniform perpendicular U{2) magnetic 
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field. We present the spectra for two different geometries: tfie plane and tlie spfiere. It 
turns out tliat tliis Hamiltonian can be mapped exactly to the one of a 2D electron in a 
perpendicular U{1) magnetic field, when the Thomas term is present and an additional 
perpendicular U{1) electric field E is applied. 



3.1. On the plane 



In order to set this problem on the plane, we consider a perpendicular, uniform magnetic 
field B = B^Uzi where 5^ is a 2 x 2 hermitian matrix. We can always choose a basis 
where the matrix is diagonal: 

2 



B. 
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We introduced the pure number /3 = involving the magnetic length Im = 
B > 0). From now on we work with h = 1. This magnetic field is a U{2) matrix, and 
is a superposition of a U{1) field B and a SU{2) field 2f3^Baz. In view of the previous 
discussion, there is an ambiguity in the notion of a non-Abelian uniform magnetic field, 
and one has to specify the non-Abelian part of the potential A. The first kind of 
potential A = ^B x r boils down to an Abelian U{1) x f/(l) gauge group, and the 
physics is simply that of two non-interacting species of particles coupled to different 
Abelian magnetic fields. The second kind however, a constant and non-commutative 
potential, is much more interesting and leads to new physics [HI [T71 1^ : 
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(12) 



The Hamiltonian describing a particle confined to a plane in this non-Abelian 
background 

1 



2m 



(p-AY 



(13) 



enjoys the translation symmetry of the plane. Since the non-Abelian part of A is 
uniform, the magnetic translation operators are insensitive to the non-Abelian part, 
and have the usual expressions: 



-idx 



y 

2/™ 



T 



-idy + 



X 

2P 



which implies immediately the Abelian Aharanov-Bohm effect: 

,{a X h) ■ Uz 



[a-T,b-T] 



-I . 



(14) 



(15) 



The r.h.s. is simply the flux of the Abelian part of the magnetic field Bl through the 
parallelogram delimited by the vectors a and 6, and the (Abelian) magnetic length scale 
is Im- Only the Abelian part of the magnetic field is quantized, and the number of states 
in a given Landau level will only depend on the Abelian field strength B. 
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It turns out that the problem of a particle in such a non-Abelian external field can 
be mapped exactly to the Hamiltonian of two-dimensional electron in the presence of 
both Rashba and Dresselhaus spin-orbit interaction [21], and it has been first solved in 
this context by Zhang [Mj . 

Having in mind to solve this problem on the sphere in the next section, we demand 
rotational symmetry around the z axis, and we focus on the symmetric gauge: 
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This symmetric case correspond to the absence of Dresselhaus interaction in 



(16) 



and 



the Hamiltonian (13) is much simpler to solve in this case. Moreover it can be mapped 
to a Thomas term ^ in the presence of a perpendicular uniform electric field E oc [5'uz ■ 
This Hamiltonian can be expanded as: 

1 



H 



UJr. 



a) a 



V^/3(aV_|_ + acr_ 



+ 2(1 + 2/3'; 



(17) 



where a, are the usual annihilation and creation operators appearing in the Landau 
problem (see Appendix A). Up to a change of spin basis U = this is nothing but 



the celebrated Jaynes-Cummings Hamiltonian, and it is straightforward to obtain its 
spectrum: 

/I 



En 



Et 



n ± 



2/3% + ^ + 



(18) 
(19) 



3.2. On the sphere 

It can be rather instructive to solve such a problem on a sphere instead of the plane. 
The surface of the sphere being finite, the degeneracy of the Landau levels becomes 
finite too, which is very interesting for numerics. Moreover the translation invariance 
of the plane is promoted to the rotational symmetry of the sphere, and the spectrum 
decomposes into SU(2) multiplets. In the Abelian case this geometry was first solved in 
[35], and later used by Haldane [26] in the context of the quantum Hall effect. In order 



to fix our notations we recall these main results in Appendix B 



3.2.1. Field configuration A uniform perpendicular magnetic field implies the presence 



of a magnetic monopole at the center of the sphere. In the Abelian case (see Appendix 
b|, the corresponding potential A^f, must have a singularity (Dirac string) somewhere 
on the sphere, for instance at the south pole 6 = n: 

1 - cos{e) 



A 



Ab 



2 rsin(6') 



(20) 
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When Dirac's quantization condition G Z is satisfied, this singularity has no physical 
consequence as it can be moved around through gauge transformations |36] . This is the 
well known quantization of the magnetic fiux piercing the sphere, which must be an 
integer number of fiux quanta: 

B-dS = 2'KN^ i.e. Br = ^. (21) 

2r^ 



To this U{1) potential (20) we add the following SU{2) component: 



A{a) = Am> + J-A^ . (22) 



Once again this correspond to a Thomas term ([o]) with a radial uniform electric field 
E oc awr, and this will turn out to be the correct extension of the symmetric gauge on 



Br = ^,-2a{l-a)'—^ . (23) 



the plane (16). The corresponding magnetic field is 

It is quite remarkable that the strength a of the non-Abelian field need not be quantized, 
as the vector potential has no singularity on the sphere: a can be any real number. 
As the radial f/(l) field is created by a magnetic monopole, it is not very surprising 
that the SU (2) counterpart involves a non-Abelian monopole. Indeed, the potential 
^ = is nothing but the large distance asymptotic of a true non-Abelian monopole 
[25] . and in this context it is well known that there is no need for a Dirac string, nor is 
there a singularity anywhere on the sphere. 

3.2.2. Hamiltonian and spectrum The details about the derivation of the spectrum of 
a particle confined to a sphere of radius r in this non-Abelian background can be found 
in 



Appendix C The Hamiltonian is the following 



Hia) = ^ [rx (p- Aia))]' , (24) 
/mr^ L V /J 

which is a scalar under global rotations generated by J = L + S", where L is the angular 
momentum and S = is the spin. Its eigenstates form SU{2) multiplets corresponding 
to the decomposition of the Hilbert space into irreducible representations of J 

^ = (jo)©2(ji)©2(j2)©---©2(j„)©--- , (25) 

where j„ = ^"^^^ + n. The corresponding eigenvalues are (see figure [l]): 

EM = ^ - 2«(1 - «) ) , (26) 




E^{a) = ^{N^ + n)- 2a{l - a) ± ^ {2a - if n {N^ + n) + j . (27) 

As can be seen in figure [l| multiple level crossings are observed. This also occurs in the 
planar case ^3], which can be recovered from the sphere in the limit of infinite radius. 
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Figure 1. Band structure for the non-Abelian Landau problem on the sphere as a 
function of the non-Abelian field strength a. This is the case of N,^ — 7 Abelian flux 
quanta, and only the lowest part of the spectrum is shown. 



3.2.3. Recovering the plane The non-Abelian A field we considered on the sphere is 
indeed the correct extension of the planar symmetric gauge (16). The planar problem 
is recovered by taking the sphere radius r — t- oo while keeping constant the gauge field 
strength on the surface: 



2r2 



5, 



« 

— ~ p = — = 
r Jb 



(28) 



The vector potential and the magnetic field become in this limit: 
B 
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B ^ B {1 + 2/3V^) . (29) 



It is straightforward to check that the eigenvalues of the Hamiltonian E^{a) behave in 
the planar limit as: 

'1 



(30) 
(31) 



reproducing the planar spectrum. Moreover one can expand the Hamiltonian on the 
sphere in terms of L 

2 



H 



2mr^ 



+ 2a{ L 



N.},r 



and using the Holstein-Primakoff representation: 
L+ = 6V^<i> + 2ata - 6t6 , 



(32) 



(33) 
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— a) a 



ViV$ + 2ata 
2 

we recover the planar Hamiltonian: 

1 



H Ur 



a"^a + - + V2(3 [aa + aV+) + 



(34) 
(35) 

(36) 



This is not surprising in view of the mapping of this problem with a spin- 1/2 
electron under an effective non-Abelian potential coming from the Thomas term with 

— * 

a perpendicular electric field. Indeed, the infinite radius limit of a sphere in a radial E 
and B field is clearly the plane under perpendicular E and B fields. 

In summary, we analyzed the non-Abelian Landau problem on the sphere and 
obtained the spectrum. The degeneracy of Abelian Landau levels is preserved, and the 
number of states per area remains 1//^, as we expected from the absence of a non- 
Abelian Aharanov-Bohm effect on the plane. Indeed, the Abelian Aharanov-Bohm (15) 
effect enforces a minimum surface of £ for the wavef unctions. 



4. Adiabatic insertion of non-Abelian flux 



For the non-Abelian gauge group SU (2), there is no notion of a 'quantum of non-Abelian 
flux'. To appreciate this fact, we recall that in an Abelian gauge theory the number of 
flux quanta is nothing but the winding number of the potential around the flux tube, 
and this integer cannot be continuously deformed. The homotopy group of U{1) is the 
set of integers Z. However, the situation is different for the gauge group SU{2), which 
has trivial homotopy group so that any flux insertion can be smoothly deformed away. 
This explains that there is no non-Abelian analog of the notion of flux quantum, in 
agreement with the absence of quantization of the non-Abelian field strength on the 
sphere. 

Nonetheless, we wish to consider the insertion of non-Abelian flux in a quantum 
Hall (qH) fluid. In the Abelian case the celebrated Laughlin argument shows that the 
insertion of a quantum of Abelian flux leads to the accumulation of electric charge ±z/e, 
with u the filling fraction of the qH liquid. Here we wish to start from an (integer) qH 
fluid and to insert, at given location O, a non-Abelian field configuration SA centered 
at O and chosen such that (i) it generates no magnetic field away from O and (ii) 
it can be removed by a gauge transformation. After this process we perform a gauge 
transformation on the evolved state, in such a way that it lives in the same Hilbert space 
as the state we started off with. 

Before entering into details we briefly sketch the well known argument for the 
Abelian discussed by Laughlin [27] and Halperin [37]. This leans heavily upon 

the existence of a gap separating the ground state of the qH liquid from its excited states. 
Take a qH droplet and adiabatically insert a flux quantum at the origin, SA^ = This 
value of the flux is very special in the sense that it creates no Aharanov-Bohm effect and 
therefore it can be gauged away. This means that after this adiabatic insertion one can 
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map the system back to its original field configuration. The net effect of this process is 
a charge transport away from the origin leaving behind a quasihole, as single particle 



states in the lowest Landau level (see Appendix A) evolve according to 



|m) -> |m + 1) . (37) 

The usual description of adiabatic processes involves Berry phases [3HI [33 , or more 
generally Berry matrices |10] in the presence of ground state degeneracies. Before we 
specify our system in detail we first give a quick review on the derivation of the Berry 
matrix. 

4.I. Berry matrix 

Considering an adiabatic process in the presence of energy degeneracies the Berry phase 
should be generalized as was done in [IQ] . We will give the details on the Berry matrix 
that are needed for our purposes. Suppose that the Hamiltonian is a smooth function of 
a parameter {\{t)} and at every point in parameter space has a degenerate ground state 
energy separated from higher levels by some finite gap. Starting out in a state |a(0)) 
belonging to the lowest energy subspace of the total Hilbert space 'He{q), the adiabatic 
theorem tells us we end up in an eigenstate which is again an element of the subspace of 
ground states at time t, \a{t)). During this process the eigenstates obey the Schrodinger 
equation H{\)\a{X)) = Ea{X)\(y{X)) . Since nothing forbids this state to pick up a phase 
or a unitary matrix the final state can be written as 

\Mt)) = e-t/o^oWf/s(t)|a(t)), (38) 
where Ub is the Berry matrix generalizing the Berry phase. It is a unitary mapping 
UBit) : T-LEoit) — > '^Eoit) from the subspace of ground states to itself such that Ub{0) = I. 
Such a map is known as a holonomy, cf . [39] . The phase in front of the Berry matrix is 
the dynamical phase, depending on the evolution of the ground state energy, but it will 
not be of any importance for our consideration, so we will discard it. The Berry matrix 
can be written in terms of path ordered integrals 



UB{t) = Vexp 



i [ A{t')dt' 
. Jo 



(39) 

l + ^e' dtn dtn-l... dtiA{ti)...A{tn), (40) 

Jo Jo Jo 



where Aa,i3 = i{a(t)\^\f3(t)) is the Berry connection, which behaves under unitary 
transformations as a gauge potential. 

4-2. Choosing the field configuration 

We consider a non-relativistic spin-1/2 particle confined to the plane in a perpendicular 
magnetic field, = BI. Writing the vector potential we use cylindrical coordinates 
and choose the symmetric gauge 
Br 

= ^I- (41) 
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The Hamiltonian of the system H = ^{p — ^'^4)^1 acts identically on both spin states, 
i.e. the Landau levels are doubly degenerate. We only consider the lowest Landau level 
(LLL) and write for the eigenstates \m,e) where e G {t?^}- For more details on these 
eigenstates we refer to Appendix A To derive the effect of inserting a field configuration 
6A{r) in this system, we will consider two particular field configurations labeled as 
M = and M = —1, respectively. In Appendix D| we present a generic configuration. 



labeled by an integer M, of which these two are specific cases. There we also give a more 
detailed derivation, to avoid any cumbersome equations in the main body of this paper 
and we explain how the label M can be interpreted. Mimicking the insertion of Abelian 
flux briefly mentioned at the start of this section, we will insert a gauge field in such a 
way that no additional magnetic field is created away from the origin. Furthermore, we 
choose a symmetric gauge and make the simplification dz{SA) = 0. The field now looks 
like a pure gauge 

6A = iU{X)VU\X), (42) 

for some unitary matrix U{X), which depends on a parameter A controlling the adiabatic 
process. The evolved Hamiltonian is now easily found to be 

H{X) = U{X)H{0)U\X). (43) 

Since this is just a gauge transformation we automatically meet the requirement that 
there is a gap of hwc separating the subspace of ground states from excited states at 



every point in parameter space. Also the eigenstates of (43) are easily found to be 
U{X)\m, e). This is all the information we need to construct the Berry matrix and find 
the evolved state. We separately present the results for the gauge configurations with 
M = and M = -1. 



4-2.1. The case M = 0. Our M = non- Abelian field configuration reads as follows 

(44) 



SAW = TTit)^^4 



SAW - iT(W^^ + i+iW^'' 
where we introduced cylindrical coordinates and ar = B-Ur and = a-u^. The reason 



why we label this field by M = is stated in Appendix D and will become especially 
clear for the case M = — 1. Note that for r ^ 1/A this field configuration does not 
depend on A and behaves as 

M<^~--a, + 0(l/r2) . (45) 
r 

The structure of this field corresponds to shifting the orbital of a spin-'|" (spin- J,) particle 
by +1 (—1), precisely what would happen by inserting an Abelian flux quantum, where 
the sign depends on the spin of the particle. From this point onwards, we will refer to 
such a field as a a^-fiux quantum, to explicitly distinguish between it and an insertion 
of a purely Abelian flux. 

Our starting point is a fully polarized integer qH state, represented as a product 
state |'?/'(0)) = (S)mLo 1"^ t); which has the the first (mj + 1) LLL orbitals filled with 
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Figure 2. Mixing coefficients Mm (A) and Um(A), for two different values of A. Around 
m ^ 1/(2A^), Um,(A) and w„i(A) are equal to each other. 



spin- 'I" particles. We adiabatically insert the non-Abelian flux (44) by slowly sweeping 
A from A = to its final value. Just like in the Abelian case we gauge this evolved state 
back to the initial situation. In this way the final state lives in the same Hilbert space 
as the initial one, 1-Le{q)- The resulting final state is as follows 

iV^o(A)) = (g) {u^{\)\m t) - Vrn{X)\m + U)) . (46) 

m=0 

The mixing coefficients {um(A), fm(A)} depend on both the orbital m and the adiabatic 



parameter A. They are given in (D.8) and are plotted in figure [2j as a function of m, for 
two values of A. Around orbital number m ~ 1/(2A^), Um{\) and fm(A) cross, resulting 
in a vanishing of the 2;-component of the spin. 

The asymptotic behavior of the mixing coefficients as m — t- oo can be read off from 



(D.9). For large m, corresponding to a radius r ^ 1/A, the adiabatic process boils down 
to |m t) — > \m + 1 \). This follows directly from the Berry matrix calculation, but it 
can be understood simply from the conservation of = + ^ctz. At large distance 



the flux we insert is essentially a a^-fiux quantum ( 45 ) , and it induces a charge transfer 
\m) — )■ |m + 1), since all particles are spin- 'I", changing the angular momentum Lz by 
one unit. Then the only way to accommodate the conservation of is through an 
accompanying spin flip | f) — )• | |). 

We can analyze the effect of this adiabatic insertion on the product state, by looking 



at the density and spin profile of the final state (46). The density is given by 



p(r;A) = y — [umiX? + v^{\f- r , (47) 

m=0 ^ \ ' J/ 

and is shown in figure |3] for four different values of A. The solid line shows the droplet 
before insertion of the non-Abelian field configuration, this is a fiat profile. Once we 
insert flux, charge is depleted from the origin and deposited at the edge of the droplet. 
This is exactly one unit of charge. 
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Figure 3. Density profile for different values of A of a product state where the first 
150 orbitals are filled with spin-^ particles. Before flux insertion (A = 0), there is a 
flat profile, but for finite values of A, a quasihole of unit charge is created around the 
origin. 



The expectation value of spin in the 2;-direction of these configurations is depicted 
in figure |4j The state before flux insertion is the blue solid line, which has trivial spin- 
texture. Upon increasing A, particles move one orbital out while flipping their spin. 
This motion starts at the outer edge of the sample and propagates towards the center. 
When the final value of A has been reached, the particles constitute a spin-texture of 
size 1/A with spin-'|" at the origin and spin-J, at the edge of the droplet. Figure [s] shows 
the spin field after inserting a flux parameterized by A = 1/3. The right picture displays 
the (x, ?/)-components of the spin field, showing that the spins have in plane winding 
number 1. 

The charged spin-texture created by the insertion of non-Abelian flux is recognized 
as a quantum Hall Skyrmion of unit electric charge, q = e, and unit topological charge, 
Qtop = 1- The topological charge, given by the Pontryagin index Qtop, measures the 
winding of the spin vector around the system (see, for example, [H] chapter 7) and can 
easily be determined by looking at figure [5j The left figure shows that the spin in the 
z-direction points up in the origin and down at the edge, so the in plane winding cannot 
be deformed into a trivial texture. 



4-2.2. The case M = —1. Our second non-Abelian field configuration corresponds 
to the M = — 1 case of the generic fiux presented in [Appendix D We will insert it 
adiabatically into the initial setting (41). The gauge potential reads 

6Ar{X) 



(4J 



There is a subtlety which did not arise in the previously discussed configuration and 
which will shed light on why we label the different fields by an integer M. This time 
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Figure 4. Expectation value of the the z-component of the spin field in a product 
state where the initial state has its first f 50 orbitals filled with spin-f particles. These 
profiles are after insertion of non-Abelian flux, for different values of A. An insertion of 
finite A creates a non-trivial spin-texture. The radius at which the expectation value 
of (Tz equals zero is around r ^ 1/A. 
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Figure 5. Spin-texture obtained after inserting non-Abelian flux with A = 1/3. The 
three components of the spin field are shown in the left picture. In the origin the spin 
points up and at the edge it points down. The figure at the right depicts the x- and 
y- components of the spin field. These two figures clearly show that a Skyrmion with 
in-planc winding number 1 is created. 



SA^{0) = -yOz 7^ 0, which is the insertion of a cr^-flux quantum, resulting in a shift of 
orbital number depending on the spin of the particle 

|m t) -f 1 t) , |m I) |m - 1 |). (49) 

The adiabatic process consists of two parts now. We start by adiabatically inserting 
a Gz — flux quantum, leading to the configuration (48) at A = 0. After that we slowly 
sweep A so as to reach its final value. Note that at every point of the adiabatic process 
we are able to find the eigenstates of the evolved Hamiltonian. Again starting from a 
product state of spin-f particles and gauging back to the original Hamiltonian after the 
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Figure 6. Mixing coefRcients of ( 50 ) for two different values of the adiabatic parameter 
A. The point where Um(A) = Wm(A) is around m ^ The asymptotes of the 

coefRcients are exactly opposite from the M — Q case. 



adiabatic process we get a final state 



l^-i(A)) 



(Mm+i(A)|m + 1 1) - Wm+i(A)|m i)) 



(50) 



m=0 



where the coefficients can be found in (D.8) and are plotted in figure |6| This time the 
scale at which the spins are ffipped is set by r ~ A. 



The density of (50) is shown in figure [7| before flux insertion and for the values 
A = 0, 1, 10. Again charge is depleted from the origin, but this time the depth of the 
hole is largest for A = 0, meaning after the insertion of a a^-flux quantum at A = 0. 
Upon increasing A, particles move inward while flipping their spin. This motion starts 
at the origin and moves out towards the edge of the sample. The resulting spin-texture 
is depicted in figure [s] and figure [oj The electric charge (g = e) and topological charge 
Qtop = 1 agree with those of the Skyrmion found for M = 0. 



5. Conclusion 

In this paper we analyzed and solved the non-Abelian Landau problem on the sphere 
and analyzed the charge and spin dynamics induced by the insertion of non-Abelian 
flux in an otherwise Abelian background. 

We remark that, in the usual (Abelian) quantum Hall setting, quantum Hall 
Skyrmions arise due to a balance between the effects of the Zeeman energy, which favors 
single overturned spins, and the Coulomb interaction, which favors configurations with 
small spin-gradients [l2jj. It is quite remarkable that our procedure of driving the non- 
interacting polarized electron gas with non-Abelian external flux leads to the very same 
Skyrmion configurations. 

Repeating the non-Abelian flux insertion in a background of a = 2 integer 
quantum Hall state, with both the spin-'l' and spin- J, lowest Landau levels completely 
ffiled, has a very different effect. In this case, the bulk state cannot accommodate 
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Figure 7. Density of a product state before flux insertion and after for different values 
of A. The flat profile depicts the initial product state. For A = a (T^-flux quantum is 
inserted creating a quasihole at the origin. At finite A, there is still a density depletion 
around the origin, but it is less sharp. 
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Figure 8. Expectation value of the z-component of the spin field, in a product state 



(50) for different values of A. The radius at which the state is unpolarized is around 
r~ A. 
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Figure 9. The left figure shows the spin field for a final state labeled by A = 5. In the 
origin the spin is pointing down, at the edge it points up. The right figure shows x- 
and y-components of the spin field, from which we see that this flux insertion created 
a Skyrmion with in-plane winding number —1. 
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any spin-flips and the effects of the flux-insertion are limited to the edges. Inserting 
non-Abelian flux through the central hole in a Corbino disc leads to neutral 5*2 = ±1 
excitations at both inner and outer edge. This situation is in many ways reminiscent 
of a thought experiment where a minimal amount of Abelian flux inserted into a 2D 
spin quantum Hall (SQH) topological phase acts as a spin pump, resulting in neutral 
Sz = ±1/2 excitations at the edges |28] . 

The details of the charge and spin dynamics associated to the insertion of non- 
Abelian flux depend on the specific form of our gauge potentials and on the way these 
depend on the sweep-parameter A. One expect that many features, in particular the 
topological quantum numbers characterizing the resulting spin-textures, will be robust 
against changes in the detailed shape of the external gauge potentials. 



Appendix A. Landau levels on the plane 



For the purpose of being self-contained, and also in order to fix notations, we recall 
the main results about the Landau problem on the plane. We consider a particle of 
charge q and mass m confined to a plane, under an external perpendicular magnetic 
field B = Buz (with qB > 0). The standard choice for the vector potential A is: 



A 



B 
2" 



/-y\ 



X 





Br 



(A.i; 



This is called the symmetric gauge because it behaves as a vector under rotations around 
Uz- The only scale of the classical problem is the cyclotron frequency Uc- 

qB 

m 

The quantum mechanical problem has an additional scale, the magnetic length 



(A.2) 




The Hamiltonian reads in the symmetric gauge: 

2 



H 



:UJr. 
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+ 




(A.3) 



(A.4) 



It is very convenient to go to complex coordinates (rescaled by the magnetic length), 
and to introduce two commuting families of creation and annihilation operators: 



V2(d + 

V2{d + 



at = i^-d + 
6t = (-3 + 



(A.5) 
(A.6) 



iz and 8=4=- In these notations the Hamiltonian and angular momentum 



where 5 — ^ <v — «- 

az az 

have a very simple expression: 



H 



a) a + 



a) a 



(A.7) 
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Figure Al. Density profile of few orbitals in the lowest Landau level, namely 
m = 2, 12 and 30 LLL orbitals. The cylindrical symmetry is clearly visible. 

from which the spectrum = ujc{n + 1/2) follows immediately, b and 6^ are also the 
generators of magnetic translation. Since they commute with the Hamiltonian, all these 
eigenvalues are infinitely degenerate. The subspace of energy En = ijJc{n+ 1/2) is called 
the n*'^ Landau level (LL). 

Denoting by n and m the eigenvalues of a^a and b^b respectively, the Hilbert space 
is spanned by the states \n,m) for n,m > 0. The additional quantum number m is 
related the value of the angular momentum Lz\n, m) = {m — n)\n, m). 

The explicit form of their wave functions is known and involves a special class of 
functions called Hermite polynomials. In this appendix we focus on the lowest Landau 
level n = 0: it is obtained by acting with b"^ on the state |0, 0): 

ffet)'" 1 
10, m) = ^-^10,0) (z|0,n) = ^= ^ exp(-z^/4) . (A.8) 
Vm! V27r V2'"m! 



The orbital |m) has angular momentum m, and the support of the wave function ( |A.8 ) 



is a ring located at distance v2m (in magnetic length scale) from the origin, as can be 



seen from figure Al 



Appendix B. Landau levels on the sphere 

In this section of the Appendix we set g = = 1. 

Appendix B.l. Field configuration: magnetic monopole 

On the sphere, a uniform perpendicular magnetic field B = implies the presence 

of a magnetic monopole at the center of the sphere, and the potential A must have a 
singularity (Dirac string) somewhere on the sphere. The gauge where the singularity 
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lies at the south pole, i.e. 

1 - cos(^) ^ 

and the gauge where the singularity lies at the north pole: 

- Nq,l + cos(9) ^ , 
2 rsin(^) ^ ^ ^ 

are related by the following unitary transformation U = e^^'i''^. This operator only makes 
sense on the sphere when Nq, is an integer. This is Dirac's quantization condition [36]. 

Appendix B.2. Hamiltonian and spectrum 

The Hamiltonian of a particle confined to the sphere of radius r in the background of 
such a magnetic monopole is: 

H= with ^ = r^[p-A]. (B.3) 

2mr"' V / 

The operators have the following (gauge invariant) commutation relations 

[A„, Afe] = ieabc (K + (r ■ B)r^ (B.4) 



and the generators of (magnetic) rotations have the form 

L = A-{r-B)f=L-^Ur. (B.5) 
They generate a SU(2) symmetry: 

[La, Lb] = ieabcLc (B.6) 
and the Hamiltonian boils down to the Casimir L^. Indeed the relation: 

A^ = L^-(^y (B.7) 

• ensures that all La commute with the Hamiltonian. 

• gives the spectrum of the Hamiltonian: ^/(/ + 1) — (^)' 

The last statement simply comes from the SU(2) algebra obeyed by the L^'s, which forces 
the eigenvalues of L^ to be of the form /(/ + 1) where / G |N. However not all these 
values of / are possible. Using the explicit expression of L, Wu and Yang [35] obtained 
the following decomposition of the Hilbert space into irreducible representations of the 
SU{2) algebra generated by L 

+ i)®-e(f + n)e... (B.8) 



and the (Abelian) spectrum on the sphere finally reads: 

_L_(„(iv, + „ + i, + ^ 



= ^2 ("(A'* + " + 1) + ^ ) " > . (B.9) 
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Appendix C. More details about the non-Abelian field on the sphere 



In this section we derive the spectrum of the Hamiltonian 

1 r / ^ m2 



describing a particle confined to a sphere of radius r in the non-Abelian background 
potential 

r X a 



A{a) = AAb + a- 



(C.2) 



where AAb stands for the U (1) potential (B.l ). Note that there is a gauge transformation 



mapping a — )■ 1 — a implemented by the unitary transformation U 



(Jr. 



Appendix C.l. Rotational symmetry and decomposition of the Hilbert space 

There are two sets of SU(2) generators in this problem: 

• the usual (Abelian) action on the coo rdina tes implemented by 
L = fx (p-AAb 



defined in (|B.5|) 



;cr. 



• the rotations in spin space generated by S — 

The Hamiltonian we are considering is not invariant under L and S separately. However, 
it is a scalar under global rotations generated by J = L + S", as can be seen from the 
expansion in terms of J: 



H(a) 



1 



2mr^ 



J2 + - - 2a(l -a) + (2a - 1) J ■ a 



-U 



-U 



. (C.3) 



1 

4 ' y - 2 ' 2 " J ' 2 

Therefore this Hamiltonian is block diagonal with respect to the decomposition of the 
Hilbert space into irreducible representations of J. This decomposition follows directly 



from the Abelian one (B.8): 
Appendix C.2. Spectrum 



A$ + 1 



A$ + 5 



(C.4) 



Working in the subspace = j{j + 1), we simply need to diagonalize the term 
X = {2a - 1) (j ■ a - ^ + ^u) + ^U. We first derive the following two relations 



J 



a 



(C.5) 



(C.6) 



2J 4 
The first one is a consequence of the gauge equivalence UH{a)U = H{1 — a), and the 
second one can be checked using the explicit form of J. From this we deduce that 
is a constant: 



(2a - 1)M + ^ ) + a(l - a)N, 



4 



(C.7) 
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± 



\ 



(2a- 1)^ 



J + 




(C.8) 



As can be seen in (25), for j > ^^^^ there are two representations of spin (j), and 
from (C.5) both eigenvalues A^.^ belong to the spectrum. However there is a unique 

— . In this irrep. U = I and the corresponding eigenvalue 
we get the following spectrum for the 



representation of spin j 
is the positive one: Rewriting j 
Hamiltonian: 



n + 



2mr2 
1 



2a(l 



(C.9) 



n{N^ + n) - 2a(l - a) ± ^ {2a - 1)^ n (iV$ + n) + {N^/2f \ . (C.IO) 



Appendix D. Generic non-Abelian field configuration 

In this appendix, we give a detailed derivation of the final state obtained after an 
adiabatic insertion of non-Abelian flux. This is done for the generic case of which two 
specific examples are discussed in section |4j The field configuration we insert is the 
following 

5Ar{\) = -{M + \) 



_^ (D.l) 

where M can be interpreted as the number of cx^-fiux quanta inserted which will 
be explained below (D.3). Inserting such a field boils down to performing a gauge 



transformation on the system 



1 




(D.2) 



So at every point of the adiabatic process we know the LLL eigenstates of the evolved 
Hamiltonian, they are given by 

\a{\)) = UM{X)\m,e). (D.3) 

Before we proceed with calculating the Berry matrix an important subtlety needs to be 



considered. We wish to insert this field configuration into a background (41). But at 
A = and for M 7^ (D.l ) is given by 5^4^(0) = {M/r)az 7^ 0, which means we have to 



start by adiabatically inserting M a^-fiux quanta, resulting into a shift of the orbitals 
depending on the spin of the particle 

|mt)^|m-Mt), \m I) \m + M I). (DA) 
After the insertion of these cr^-flux quanta, we slowly sweep A from zero to some final 



value resulting into (D.l). Now we can use the eigenstates (D.3) to compute the Berry 
connection 

= ^{a{X)\J^\f3{X)) = ^(a(0)|f/l^f7M|/3(0)), (D.5) 



Non-Ahelian gauge fields 



23 



where 

i\ f _^2M+l \ 

The Berry connection only has non-zero elements between states of the form 
{f^A/(A)|m t), f/M(A)|m + 2M + 1 I)}. Written in this basis, for every m the Berry 
matrix is a 2 x 2 matrix 

f^B = cos(eL^)(A))I + tsm{ei^\X))ay, (D.7) 

where the angle is given by 

roo 2+2M+2m 

el^'\X)= / rfrarctan(Ari+2M) (D.8) 

Jo 2-+A^v^2m!(m + 2M + l)! 

This angle has interesting asymptotes in two different limits 

lim el^'^ = arctan(A(2m)^^+i/2)^ ^ 

lim = - r(,n + M + 3/2) ^ 
2 ^m!(m + 2M+ 1)! 

After the adiabatic insertion of flux we gauge the system back to the initial one. This 
cycle has the following effect on a single particle state |m '|") 

Ul{X)U^{X)UM{X)\m t) = ul^'\X)\m t) - vl^'\X)\m + 2M+1 1), (D.ll) 

where the mixing coefficients are expressed in terms of ( D.8[ ) 



ul^'\X) ^ cos(eHA)) , vl^'\X) ^ sin(eHA)) • (D.12) 
That the equality in (D.ll ) holds can be seen by inserting unity Um{X) |m', e) (m', e\Ulj{X) 



between Ul,^ and U^. After deducing the effect of the two stages of the adiabatic pro- 
cess, we can combine them to find the final state. Before we give the final state, one last 
remark needs to be made. Since we want to stay in the LLL, we have to put the state 
on a Corbino disc, meaning that we fill the orbitals of the initial product state with 
spin- 'I" particles starting from some initial orbital rrii up to a final orbital ruf. The two 
specific adiabatic flux insertions given in section |4] are actually the only two situations 
for which the Corbino disc is not a necessary geometry for staying in the LLL. Starting 
from a product state on a Corbino disc where the orbitals are filled with spin-"|" particles, 
the final state after first adiabatically inserting M cr^-fiux quanta, then cranking up the 



value of A in (D.l ), and finally gauging back to the initial configuration, is given by 

IV^m(A))= (g) [ul'f\mrn-Mt)-vZ\,W\rn + l + Mi)).{D.13) 

m=mi 
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